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Abstract 

Central to many discussion of decoherence is a master equation for the reduced density matrix of 
a massive particle experiencing scattering from its surrounding environment, such as that of Joos 
and Zeh. Such master equations enjoy a close relationship with spontaneous localization models, 
like the GRW model. This aim of this paper is to present two derivations of the master equation. 
The first derivation is a pedagogical model designed to illustrate the origins of the master equation 
as simply as possible, focusing on physical principles and without the complications of S'-matrix 
theory. This derivation may serve as a useful tutorial example for students attempting to learn 
this subject area. The second is the opposite: a very general derivation using non-relativistic 
many body field theory. It reduces to the equation of the type given by Joos and Zeh in the 
one-particle sector, but correcting certain numerical factors which have recently become significant 
in connection with experimental tests of decoherence. This master equation also emphasizes the 
role of local number density as the "preferred basis" for decoherence in this model. 

PACS numbers: 03.65.-w, 03.65.Yz, 03.65.Ta, 05.70.Ln 
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I. INTRODUCTION 



Non-unitary master equations for a density matrix arise in both continuous state localiza- 
tion models, such as GRW theory [1], and in decoherence calculations in standard quantum 
theory. They differ, however, in their underlying physical pictures. GRW theory involves 
an explicitly modified dynamics in the Schrodinger equation. Standard decoherence calcu- 
lations, by contrast, employ the usual Schrodinger dynamics for a system coupled to its 
environment, but this becomes a non- unitary dynamics for the reduced density matrix once 
the environment is traced out. These similarities and differences have undoubtedly been 
clear ever since the appearance of GRW theory twenty years ago, but there is surely still 
more to learn in this area. It is a great pleasure to have the opportunity to contribute to 
this volume in honour of the 70th birthday of GianCarlo Ghirardi. 

This paper concerns master equations in standard quantum mechanics of the type used 
in decoherence studies, obtained by coupling a point particle by some kind of environment. 
A simple example of such a master equation in one dimension is 



This describes quantum Brownian motion for a free particle of mass M in the limit of 
negligible dissipation. This is the one- dimensional version of an equation first obtained 
by Joos and Zeh [2] (in the limit of small \x — y\), which involves a massive free particle 
undergoing scatterings by an environment of much lighter particles. Many other derivations 
of this and similar equations have since been given [3-7]. This equation can also arise for 
the case of a point particle linearly coupled to a thermal bath of harmonic oscillators, and 
this model also has been the subject of many papers [8-10]. Models involving fields have 
also been considered [11-13]. (The literature on quantum Brownian motion is considerable 
so only a selection is mentioned here. See also Ref. [14].) 

A key reference point in these studies is the Lindblad form of the master equation [15] , 
which is the most general possible form a master equation can take under the assumption 
that the evolution is Markovian (a condition well-satisfied in a wide variety of interesting 
models). The Lindblad master equation is 




(1.1) 
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Here, H is the Hamiltonian of the distinguished subsystem (sometimes modified by terms 
depending on the Lj) and the n operators Lj model the effects of the environment. For 
example, the master equation of one-dimensional quantum Brownian motion, including dis- 
sipation, is of the Lindblad form with a single Lindblad operator 



as described in Refs.[16, 17]. (This reproduces Eq.(l.l) for small 7 with D = 2M^kT/h 2 ). 

Equation (1.1) describes the decoherence process in which an arbitrary initial density 
matrix becomes approximately diagonal in position on a very short timescale. This process 
is thought to be a key element in understanding how classical behaviour emerges from 
quantum theory [18-20]. Recent experiments have been able to actually observe the rate 
of the decoherence process [21], which is connected to the constant D in Eq.(l.l). It turns 
out that some of the original derivations of the master equation led to incorrect values of 
D. More recent derivations [6, 11] have corrected these errors and produce values of D 
compatible with experiments. (See also Refs.[22, 23]). 

The purpose of this paper is to present two derivations of a class of master equations of 
the form Eq.(l.l) and its generalizations. 

The first derivation, described in Section 2 and 3, is a simple pedagogical model, designed 
to illustrate the way in which the general form of Eq.(l.l) follows from some simple physical 
ideas. We therefore avoid the technical complications and sometimes non-transparent math- 
ematical assumptions involved in these models, but make no claims about making physically 
accurate predictions. (This model is similar to one presented by Joos et al. [24]). 

The second derivation, described in Sections 4, 5 and 6, is a very general derivation of 
a class of master equations and makes use of non-relativistic many body field theory. It 
reduces, in form, to an equation of the type given by Joos and Zeh and others, but with a 
correct value of the decoherence rate that is compatible with experiments. This derivation is 
also relevant to another issue in decoherence theory, which is the question of the "preferred 
basis" - the natural basis in which interferences are destroyed. In the master equation 
Eq.(l.l), it is clearly the interferences between different values of position that are initially 
suppressed, so that position is the preferred basis. More generally, it is known that under 
evolution according to the Lindblad form Eq.(1.2), it is the Lindblad operators Lj that 
defined the preferred basis (at least in simple models) [25]. We will see in the many body 




(1.3) 
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field theory model that most generally in system-environment models, local number density 
is the preferred basis (with position emerging as a special case of this in the one-particle 
sector). The special role of number density also echoes certain aspects of the GRW model 
[1]. Furthermore, local densities are thought to play a key role in the most general possible 
derivations of emergent classicality, even when there is no environment present [20, 26] 

The many body field theory derivation of the master equation first appeared, in essence, 
in Ref.fll], as part of a wider investigation into the properties of the decoherent histories 
approach to quantum theory [27-31]. Here it is presented on its own merits as a contribution 
to the theory of quantum Brownian motion. 



II. A SIMPLE MODEL 



We consider a system consisting of a particle in one dimension which interacts through 
occasional collisions with a gas of light particles. Except for the collisions, which are assumed 
to be very brief, the particles evolve freely. By considering the change in the system density 
matrix during these collisions we will derive the form of the decoherence term (the last term) 
in Eq.(l.l). 

We first consider the collision process classically. Suppose the system particle has mo- 
mentum P and mass M and a particle from the environment has momentum p and mass 
m. We will assume that the collision conserves both energy and momentum. If the final 
momenta are P' and p', we therefore have 

p2 2 p,2 ,2 

2M 2m 2M 2m K ' J 

P + p = P' + p' (2.2) 

Ignoring the trivial solution P — P' and p = p', the final momenta are given by 

(M-m) 2M 

~ (M + m) + (M + mf ( 3) 

= aP + bp (2.4) 

/ 2m (M-m) 

P = (M + m) ~ (M + mf (2 ' 5) 

= cP-ap (2.6) 

(where the coefficients a, 6, c are read off from Eqs.(2.3),(2.5)). It is now very useful to make 
two approximations. We first assume that the environment particles are much lighter than 



the system particles: 

m < M (2.7) 

Second, we assume that the speed of the system particle is much smaller than the speed of 
environment particles, 

P v 

77 < - ( 2 - 8 ) 
Mm v ' 

(although the momenta P and p may be comparable in size). These approximations imply 
that the momenta after the collision are given by the much simpler expressions 

P' « P + 2p (2.9) 
p' « -p (2.10) 

Turn now to the quantum case. We will make the reasonable assumption that energy and 
momentum are also conserved by the quantum description of the collision. The key idea in 
the quantum case is to work with states of definite momentum (plane wave states) and use 
the above results to deduce how they change during a collision. We will eventually also need 
to assume something about the locality of the interaction, but that will not be needed just 
yet. 

Introducing the positions (x, q) of the system and environment particles, it follows from 
Eqs.(2.9), (2.10) that an initial plane wave for the total system changes according to 

e iP* e fa _> e ^ p+2p)x e-^ q (2.11) 
= e * Px e^ 2x -ti (2.12) 

(in the approximations described above). Or, in configuration space, the effect of the collision 
is to make the replacement 

x -> x (2.13) 
q -> 2x -q (2.14) 

Since the plane waves are a complete set of states, this result determines the effects of a 
collision on any initial state. 

We suppose that the initial density matrix of the whole system is a simple product state, 
p(x,y) ps(q,q')- As a result of a single collision, Eqs.(2.13), (2.14) imply that the total 
density matrix changes according to 

P(x, y) pe(q, q') -> p(x, y) p £ (2x - q,2y - q') (2.15) 



Tracing over the environment, this means that the density operator of the system only 
evolves according to 

p(x,y) -> p(x,y) J dq p e (2x - q,2y - q) (2.16) 

We suppose that the collisions take place at a rate Y per unit time. The total change in the 
density matrix of the system during a small time interval At is therefore given by 

Ap(x, y) = -At F(x, y) p(x, y) (2.17) 

where 

F(x, y) = T (l- J dq p £ (2x -q,2y- qfj (2.18) 
This is also usefully written, 

F(x, y) = T Tr ((1 - e 2 ^ x ~^)p^ (2.19) 

where p is the momentum operator on the environment and the trace is over a complete set 
of environment states. It follows that the master equation is 

a P m (p P d>p\_ F{XtV)p (2 . 20) 



dt 2M \dx 2 dy 2 
In the small \x — y\ limit we have 

F(x,y) = -^±(x-y) 2 (2.21) 

where we have assumed that the environment state is such that (p) = 0. For a thermal 
environment state we have 

(p 2 ) = mkT (2.22) 

so we obtain a master equation of the form Eq.(l.l) with D = 2mTkT/h 2 . We have therefore 
obtained the expected form of the master equation, using energy and momentum conserva- 
tion, together with the approximations that the environment particles are much lighter and 
faster than the system particle. 

A slightly different but simpler model along these lines was given by Joos at al [24] . Their 
model postulates a simple dynamics involving a particle being subject to random kicks, but 
without relating it to a collision process with energy and momentum conservation as here. 
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III. A WIGNER FUNCTION DERIVATION INCLUDING DISSIPATION 



It is now useful to give a more detailed derivation of the master equation in this simple 
model using the Wigner representation. We go beyond the approximations used above and 
work to leading order in m/M. This derivation shows how the dissipative terms arises and 
also establishes the connection between the collision rate Y and the dissipation 7, hence 
connects microscopic and macroscopic parameters. 

The Wigner function of the density matrix p(x, y) of a one-dimensional systems is defined 

by 

W(p, ?) = ^ / # ^ P(Q + \*><l- \t) (3-1) 
together with its inverse 

p(x, y) = j dp W(p, ^) (3.2) 

The Wigner function has the properties 

dpW(p,q) = p(q,q) (3.3) 
dqW(p,q) = p(p,p) (3.4) 



/ 



where p is the Fourier transform of p(x, y), so W contains the usual position and momentum 
probabilities as its marginal distributions. The Wigner function is the closest thing quantum 
mechanics has to a phase space probability distribution function but narrowly fails since W 
is not always positive. Its time evolution is identical to classical evolution for linear systems, 
with corrections proportional to powers of h 2 for non-linear potentials. It is therefore a very 
useful tool for discussing the connection between quantum and classical systems, although 
note that there are subtle differences between the Wigner function and a classical distribution 
function, as we shall see shortly. (See Refs.[32] for properties of the Wigner function.) 

In the absence of interactions, the two particle Wigner function W 2 (P, X,p, q) obeys the 
equation 

aw 2 = Pdw 2 p dw 2 

dt ~ M dX m dq ^ ' ' 

and this corresponds to unitary evolution of the density matrix. We need to find terms 

representing the collision to add to the right-hand side, using an argument along the lines 

of that used in the previous section for the density matrix. We again consider the collision 

process described by Eqs.(2.4), (2.6) and look for the change in the Wigner function. 



If W were a classical distribution function, it would be sufficient to consider only the 
momentum transfer described by Eqs.(2.4), (2.6). In the Wigner function, however, things 
are a bit more subtle. If Eqs.(2.4), (2.6) describe a process in the momentum representa- 
tion of a quantum system, then there is a corresponding transformation in position space. 
Proceeding along lines identical to the derivation of Eqs.(2.13), (2.14), it is easy to see that 
this transformation is 

X -> aX + cq = X + c(q-X) (3.6) 
q -> bX - aq = 2X - q + c(q - X) (3.7) 

That is, the collision process in the Wigner function involves both the transformation 
Eqs.(2.4), (2.6) on the momenta together with the transformation Eqs.(3.6), (3.7) on the 
positions. This must be the case because on integrating out the momenta in the Wigner func- 
tion, the correct distributions for position must be obtained, and vice versa (as in Eqs.(3.3), 
(3.4)). 

Denoting the system Wigner function by W(P, X) and the environment Wigner function 
by Ws{p, q), the above discussion implies that the effect of a collision is to produce the 
transition 

W(P, X) Wsip, q) -> W(aP + bp, aX + cq) W £ (cP - ap, bX - aq) (3.8) 

in the two particle Wigner function. However, in this derivation it turns out that it is 
also important to incorporate the fact that the interaction is local, that is, it is described 
by a potential of the form V(x — q) which decays for large \x — q\. In a more complete 
derivation this would be accomplished by multiplying Eq.(3.8) by a function of X — q which 
is concentrated around X — q. However, in the interests of keeping the derivation heuristic, 
we will incorporate this "by hand", by assuming that terms of the form \q — X\ are small 
in some sense. Since the constant c is also small (of order m/M), the most minimal and 
natural implementation of locality is to simply drop the terms c(q — X) in Eqs.(3.6), (3.7). 

Integrating out the environment, and considering Y collisions per unit time, we now have 
that the evolution equation of the system Wigner function is 

dW P dW 



+ 



m MdX ' T J dpdq [ W{aP + bp > X) We{cP ~~ ap > 2X ~^- W ( P > X ) W £<P> 9)1 

(3.9) 
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which is a Boltzmann equation. Using some simple changes of variable for p and q in the 
first term, this is easily rewritten, 



dW P dW „ . , , 
1H=-M8X +T ' dpdq 



-W(— - -p,X) - W(P,X) 
a a a 



W £ (p,q) (3.10) 



(From this form we see that the position coordinates in the Wigner function are, in the end, 
effectively unchanged as a result of the collision, despite the argument above, but this is due 
to the approximation of dropping terms proportional to c(q — X)). 

We again assume that the environment is in a thermal state at temperature T. To obtain 
the more familiar form of the equation, we will assume that W is slowly varying in P so 
that it may be expanded in derivatives. We will also take m <C M and look for the leading 
order terms in m/M, so we have that 

1^=,! + ^ (3.11) 
a M-m M 

We thus obtain 

dW P dW ^[2m Tir 2m dW 1 / b\ 2 , 2x d 2 w\ , n _ 

W = -M8X +T {jf W + Jf P W + 2{a) <P> a^J (312) 

where we have dropped terms of order m 2 /M 2 and terms involving higher derivatives of W. 
Now we note that it is appropriate to identify the dissipation 7 as 

7 = ^r (3,3) 

hence we make a connection between the collision rate and the dissipation. Since b/a ~ 2 
to leading order and (p 2 ) = 2mkT, we finally obtain the result 

dW P dW d(PW) 1r ^d 2 W 

This is the expected Wigner equation for a system undergoing quantum Brownian motion 
including dissipation. As is well known, this equation describes the approach to thermal 
equilibrium. Inverting the Wigner transform, the corresponding density matrix equation is 

dp ih (d 2 p d 2 p\ , Jdp dp\ 2M 1 kT . ^ 

' ' -ihyix-y) / -/ ^—{x-yfp (3.15) 



dt 2M \dx 2 dy 2 J \dx dy J h 2 

Although often given as the master equation for quantum Brownian motion with dissipation 
[8], this equation is not in fact of the Lindblad form Eq.(1.2) and actually suffers from 
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a possible small positivity violation [33, 34]. However, it can easily be modified into the 
Lindblad form Eq.(1.2) with Lindblad operator Eq.(1.3) by addition of a term with coefficient 
proportional to 1/T, so the difference is neglible for high temperatures. 

In summary, we obtain the master equation Eq.(3.9) with dissipation using momentum 
and energy conservation to describe the collision, and taking leading and first orders terms 
in order m/M together with a simple approximation to incorporate the the locality of the 
interaction. The familiar forms Eqs.(3.14), (3.15) are obtained using the further assumption 
that the Wigner function is a slowly varying function of P (which, in the density matrix, 
corresponds to exploring the region of small \x — y\). 

IV. MANY BODY FIELD THEORY 

We now turn to the second derivation of the master equation, using many body field 
theory. This section is based on Ref. [11]. We begin by briefly reviewing the formalism 
[35, 36] We consider a set of non-relativistic system particles described by a field -^(x) 
interacting through a potential 0(x) with an environment described by a field x( x )- The 
total system is described by the Hamiltonian 



(For simplicity we set H — 1 hereafter). In this language, the number densities iV(x) and 
n(x) of the system and environment fields are 



The above relations are also more conveniently written in terms of and b^, the anni- 
hilation operators for the system and environment, respectively, and the Hamiltonian then 
is 




(4.1) 



7V(x) = V f (x)^(x) 
™( x ) = X f (x)xW 



(4.2) 
(4.3) 




(4.4) 



q 





k' 1 +k^,=ki+k 2 
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where E n = q 2 /2M, cj q = q 2 /2m, V is the spatial volume of the system (which we assume 
is in a box) and 

z/(k) = J d 3 x e~ lk ' x 0(x) (4.5) 
The Fourier transformed number densities are 

N k = ^a q a q+k (4.6) 
q 

q 

and one may see that the Hamiltonian has the more concise form 



H = J2 (£q«q«q + ^q & W + ^ Yl K^W-k (4.8) 
q k 

= H + ^ int (4.9) 



From these relations we see that the environment couples to the number density of the 
system. It is this feature of many body field theory that makes it the appropriate medium 
for the derivation of the master equation emphasizing the role of number density. 
The ^-matrix is 

S = Texp (-i J dt H int (tfj (4.10) 

where 

H int{t) = ^7 E "0OlVk(f)n-k(f) (4.11) 

and here 

Mt) = ^a q a q+k e l ^-^+ k )' (4.12) 
q 

Mt) = ^ 6 i 6 q+k ei(Wq " Wq+k) * ( 4 - 13 ) 

q 

It is enlightening to look at a simple scattering situation to determine how the environ- 
ment stores information about the system (which in turn determines the preferred basis). 
Suppose, for simplicity, that the distinguished system is classical, and consider what hap- 
pens when the environment scatters off it. Suppose the environment starts in an initial 
momentum state |k ) and scatters into a final state |k/). The scattering amplitude for this 
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process, to first order, is 

• poo ^ 

(k/l^lko) = — J Jt^^)m)^f\n- k (t)\ko) 

= — ^(k) f dt N k (t) e*^/"^)* (4.14) 
IV J 

where k = kj— k . This simple result shows that a single scattering event by the environment 
stores information about the Fourier transform (in space and time) of the number density. 
It is in this sense that the number density has a preferred status - this is the variable that is 
measured most directly by the environment. (An analagous result holds in linear oscillator 
models [37]). 

The measured variables above are of course non-local in time, involving a temporal Fourier 
transform of the number density, so cannot in fact be compatible with a Markovian master 
equation of the Lindblad form. Under a reasonable slow motion assumption, the system 
timescale is much slower than the environment timescale, and we may ignore the time- 
dependence in iV k (t), yielding 

(k f \S\k ) = ± i/(k) iV k % k/ - a*J (4.15) 

This corresponds more directly to a Markovian master equation, as we shall see. 



V. DERIVATION OF THE MASTER EQUATION 

Following the method first used by Joos and Zeh [2], we may derive the master equation for 
the reduced density operator p of the system by considering the scattering of the environment 
off the system, to second order in interactions. We assume that the system and environment 
are initally uncorrelated, so the total density operator is 

Pt = Po®P£ (5.1) 

We also assume that each scattering event takes place on a timescale which is extremely 
short compared to the timescale of system dynamics. This means that in an interval of time 
At which is long for the environment but short for the system, we may write, 

p T {t + At) = Sp T {t) S ] (5.2) 
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(where we are using the interaction picture). Expanding (5.2) to second order, the S'-matrix 
may be written, 

S = l + iU 1 -U 2 (5.3) 



where 



and 



/oo 
dt H int (t) 
-oo 



(5.4) 



U 2 = l -Jdt 1 J dt 2 T(H int (t 1 )H int (t 2 )) (5.5) 
The requirement of unitarity, S* -1 = implies that Ui = Jj\ and 



U 2 + U\ = Ul (5.6) 



We will therefore write 



where B = B\ so we now have 



U 2 = hj 2 x +iB (5.7) 



S=l + i(U 1 -B)-hjl (5.8) 



Inserting this in (5.2), we obtain 

^-At = i[U 1 - B, p T ] + U.prU, - hjlpr - l -p T U 2 x (5.9) 

We now trace Eq.(5.9) over the environment to obtain the master equation for the system 
density operator p. As is usual in this sort of model, we assume that the environment is so 
large that its state is essentially unaffected by the interaction with the system. Since the 
total density operator starts out in the factored state (5.1), this then means that, to a good 
approximation, px persists in the approximately factored form p <g> ps, and we may insert 
this in the right-hand side of Eq.(5.9) [38]. We thus obtain the preliminary form for the 
master equation 

^ At = i[Tr £ (U lP£ ) - Tr £ (Bp £ ),p] + Tr £ (u^tU^ - l -U 2 lPT - \ptU 2 ^J (5.10) 

We now work out these terms in more detail. We will employ the simple but useful slow 
motion approximation, in which we ignore the time-dependence of N-^(t). (Corrections to 
this approximation are considered in Ref. [11]). This implies that 

U x » "^$>(k)iVk XXVk 27r% q -a; q _ k ) (5.11) 



2V , 
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The important terms for decoherence are the final three terms on the right-hand side of 
(5.11). When traced, these give, 

Ti £ (u^tU, - hjlpT - \ptU 2 ^J = c(k, k') (jV k ,pjV k - ^N k N k ,p - ^pN k N k ^j 

(5.12) 

where 

c(k, k') = v(k)v(k') ^2 ^(^q - ^q-k)5(^ q ' - cj q /_ k /) (& q & q _ k & q/ fc q /_ k /)£ (5.13) 

qq' 

We will take the environment to be a thermal state, which is diagonal in the momentum 
states. It follows that 

(frqfrq-kfrq'V-k')£ « <5 q , q '-k' ^q',q-k (5.14) 

This implies k = — k', and also that the two delta- functions are the same in Eq.(5.13). We 
then interpret the square of the delta-function in the usual way, 

[5(cj q - CJ q _ k )] 2 = 5(0) 5(cJq-CJ q _ k ) 

^ 5(cJ q - CJ q _ k ) (5.15) 



We now have 



2tt 



C (k,k') = 5 k> - k 'c(k)^ (5.16) 



where 

1 

c (k) = — k(k)| 2 ^ <J(w q -w q _ k )(6t6 q _ k 6T_ k 6 q ) £ 
q 

= ^i^ k )i 2 E + ( 5 - 17 ) 

q 

The terms involving environment averages have the usual thermal form (for a bosonic envi- 
ronment), 

= e ^ q -,) _ I ( 5 - 18 ) 
where (3 = 1/kT with T temperature, and p is the chemical potential. 

The form Eq.(5.16) means that the important terms in the master equation are of the 
Lindblad form, 

Tr £ (urprU, - l -U 2 lPT - ^p T ufj = At ]Tc(k) (^N k pNl - ^N^p - ^pN^ (5.19) 
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where we have used the fact that = A_ k . The remaining two terms in Eq.(5.10) clearly 
just modify the unitary dynamics of the system. First we have 

Tr £ (U lPe ) = ^J2 K k )^k E( fo qVk>£ 2tt5K - co q _ k ) (5.20) 

k q 

Clearly from the term (fr q & q - k )£ this expression will be zero unless k = 0, and therefore it is 
proportional to N, the total particle number operator (although the overall coefficient will 
need to be regularized). This therefore contributes a term to the master equation of the 
form [N, p\. We assume that there is a fixed number of system particles so it is reasonable 
to take this term to be zero. 

The other remaining term in Eq.(5.10) involves the time ordering terms in £/ 2 and is a 
bit more complicated to evaluate. Fortunately, the detailed form of this expression is not 
needed here, and it can in fact be easily shown that this term has the form 

Tr e (Bp e ) = AtJ^ N k Nl (5.21) 

k 

for some coefficient d(k) which we will not need. Inserting all these results in Eq.(5.10), the 
factors of At all drop out, and we obtain, in the Schrodinger picture, 

d j- = -i[H - rf ( k ) N * N l P\ + E c ( k ) ( N *PNl ~ \NlN v p - \pNINtS (5.22) 
k k ^ ' 

As desired, this is the Lindblad form with the Lindblad operators given by 

L k = C 5(k)iV k (5.23) 

We have therefore produced a derivation of the master equation for a scattering environment 
which shows very clearly the key role of local number density as the preferred basis, as 
indicated by the simple scattering calculation, Eq.(4.15). 

It is interesting to note that the decoherence effect is second order in interactions, but we 
were able to anticipate it from the simple first order calculation, Eq.(4.15). The reason for 
this is the relationship Eq.(5.7), which shows that the important part of the second order 
terms is the square of the first order terms, and this is a consequence of unitarity. 

VI. COMPARISON WITH PREVIOUS WORKS 

It is now important to check that the master equation we have derived reproduces known 
results when we restrict to the one-particle sector for the system, where a number of deriva- 
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tions have been given in a quantum mechanical framework [2-6]. In the one-particle sector 
we may work with a density matrix p(k, k') = (k|p|k'), or equivalently p(x, y) in the position 
representation. We use the relations 

\A] = a U (6.i) 

[iV q ,a k ] = -a k + q (6.2) 

These relations imply that 

iV q p(k, k')AL q = p(k - q, k' - q) (6.3) 

AL q iV q p(k,k') = p(k,k') (6.4) 

p(k,k')AL q 7V q = p(k,k') (6.5) 

In the position representation, this means 

iV k p(x, y)iV_ k = e tk ( x - y V(x, y) (6.6) 

The master equation for the one-particle density operator p(x, y) is then 

dp(x, y) 
dt 

where 



= -i(x| [H , p] |y) - F(x - y)p(x, y) (6.7) 



F(X " y) = jhf j d3<ld3k |Z/(k) ' 2 nq(nq_k + 1} 6{iUn ~ (1 " (6 ' 8) 

Note that the term involving the coefficient d(k) in Eq.(5.22) drops out because [iV k iV^, p] = 
in the one-particle sector. 

Eq.(6.7) is of the same general form as earlier results [2-4]. To compare in detail, we first 
introduce the quantity 

/(k,k') = |^(k-k') (6.9) 

(which appears in the usual Born approximation to first order scattering). Then, letting 
k — > — k + q in (6.8), we get 

F(r) = ^ j d 3 qd 3 k |/(q,k)| 2 n q (n k + 1) % q - LU k ) (1 - e^" 1 ^) (6.10) 

The delta-function implies that q 2 = k 2 , and we find that 

F(r) = J dq q 3 n q (n q + 1) J <ftW|/(q, k)| 2 (l - ^>') (6.11) 
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For n q <C 1, and identifying (l/2n 2 )q 2 n q dq as the fraction of particles with momentum 
magnitude between q and q+dq, we find agreement with the careful derivation of Hornberger 
and Sipe [6], which in turn agrees with recent experiments which measure the decoherence 
rate [21]. (Hornberger and Sipe corrected erroneous numerical factors in some of the earlier 
derivations [2, 3], but the qualitative and order of magnitude predictions of these earlier 
works are correct. See also the cross-check of Adler [23]). 
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